Coulomb blockade of tunneling through a double quantum dot 
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We study the Coulomb blockade of tunneling through a double quantum dot. The temperature 
dependence of the linear conductance is strongly affected by the inter-dot tunneling. As the tunnel- 
ing grows, a crossover from temperature-independent peak conductance to a power-law suppression 
of conductance at low temperatures is predicted. This suppression is a manifestation of the Ander- 
son orthogonality catastrophe associated with the charge re-distribution between the dots, which 
accompanies the tunneling of an electron into a dot. We find analytically the shapes of the Coulomb 
blockade peaks in conductance as a function of gate voltage. 



PACS numbers: 73.20.Dx, 73.40. Gk 



I. INTRODUCTION 

Electron tunneling in a mesoscopic structure may be 
significantly affected by the charging effects. The charg- 
ing suppresses tunneling, if the charge spreading is im- 
peded by weak links, or by a special geometry of the 
structure. Such a suppression of tunneling is commonly 
referred to as the Coulomb blockade, for a review see 
Rcf. |l|. In recent experimentsa it has become possible 
to observe the Coulomb blockade in semiconductor het- 
erostructures where the geometry of the system can be 
easily modified by adjusting the voltages on special gate 
electrodes. 

A common example of the Coulomb blockade effect is a 
measurement of linear conductance between two macro- 
scopic leads weakly coupled to a quantum dot.H When an 
electron tunnels from a lead to the dot, the electrostatic 
energy of the system 



U 



(1) 



changes; here C is the capacitance of the dot, en is its 
charge, Vg is the gate voltage, and k = Cg/C is a di- 
mensionless geometrical factor which defines the gate ca- 
pacitance Cg. At low temperatures, T <C e^/2C, the 
equilibrium discrete charge of the system is determined 
by the minimum of U. Tunneling of an electron into or 
out of the dot leads to a large increase of the energy, and 
conduction through the dot is suppressed. However, at 
certain values of the gate voltage the electrostatic energy 
is degenerate. 
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and the Coulomb blockade is lifted. Therefore, the linear 
conductance shows a series of peaks at the gate voltages 
V* = (2n-|-lle/2Cg. The heights and shapes of the peaks 
can be founds using the master equation technique. 
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Here Gi^r are the conductances of the weak links connect- 
ing the dot to the leads. 

In a Jiiiraber of recent experimentalcilS and 
theoreticaBoiiS papers tunneling through two coupled 
quantum dots was explored. In particular, by using a 
double dot structure one can probe the quantuttycharge 
fluctuations more directly than in a single dot.QO'EJ Here 
we focus on the geometry of Ref. ^ shown schematically in 
Fig. ^, in which the dependence of the peak positions on 
the conductance Gq of the constriction between the two 
dots was studied. We discuss the theory of the peak posi- 
tions in Sec. ||. As Go grows and approaches 2e^ /h, the 
peaks become equidistant, and in this respect the two- 
dot system behaves as a single dot of a larger size. It is 
clear, however, that unlike a large single dot, the charge 
spreading between the two coupled dots is impeded, and 
takes a relatively long time, t ^ G /Gq ~ ?i(e^/G)~^. The 
characteristic energy related to this time delay, Ti/t, is of 
the order of the charging energy, and one can expect it 
to affect the conductance through the double-dot system 
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and cause deviations from Eq. (M). Indeed, in Sec. 
show that the slow propagation of charge between the 
dots results in a suppression of the conductance peaks. 
The specific shape and temperature dependence of a 
conductance peak provides one with information about 
quantum fluctuations of charge between the two dots. 



II. POSITIONS OF THE PEAKS IN LINEAR 
CONDUCTANCE 



To discuss the Coulomb blockade, one has to introduce 
the electrostatic eoergy of the system shown in Fig. ||. 
In the experimenta the potentials of the dots were con- 
trolled by a single gate voltage Vg. Clearly, the equilib- 
rium electrostatic energy is a function of three variables: 
the discrete charges of the two dots eAi and eN2, and the 
gate voltage Vg. It also depends on the capacitances of 
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the dots — to the gate, to the external world, and to each 
other — which introduce five parameters into the problem. 
We will use the following expression for the electrostatic 
energy 

U{Ni,N2) = Ec{Ni + N2- 2Xf 

+Ec[Ni- N2 + \{Nr+N2)-aXf, (4) 

where X is a dimensionless variable proportional to Vg. 
The effective circuit we have in mind, the exact relation 
between X and V^, and expressions for the parameters 
Ec^ Ec, A and a in terms of the capacitances of the dots 
can be found in Appendix [a[ 

We intentionally grouped the terms in Eq. (^) in such 
a way that the energy depends on the total number of 
particles in the two dots Ni + N2 and the relative charge 
Ni — N2. In this paper we assume that the coupling of 
the double-dot system to the leads is extremely weak, 
Gi,Gr <C Go, and therefore one can neglect the quan- 
tum fluctuations of Ni + N2- On the other hand, the 
inter-dot conductance Go is not necessarily small, and 
at Go ~ e'^/h the fluctuations of Ni — N2 are signifi- 
cant. The tunneling of the electron between the dots 
lowers the ground state energy of the system. Thus to 
determine the positions of the peaks in the linear con- 
ductance one should generalize Eq. (||) by replacing the 
electrostatic energy U with the ground state energy of 
the double-dot system Ej^{X) for a fixed total number 
of electrons J\f = Ni + N2. That is, the peak positions 
X* are given by 



E^{X*) = E^r+l{X*). 



(5) 



Early attempts at the calculation of the peak positions 
were based pajtiodels allowing only a few discrete states 
in each dotJld'EJ Such an approach should provide an ad- 
equate description of the system in the case |Uf|-extremely 
small quantum dots. In typical experimcntsjja however, 
the number of states in each dot is large, and a model 
with continuous spectra of electrons is more appropriate. 
The calculation of the ground-state energy for such a 
model in the limits of weak and strong coupling between 
the dots can be found using the techniques developed in 
Refs. [l^ , |l5| . For a symmetric system, A = a = 0, in the 
weak tunnchng case, Go ^ e^/h, the peaks_am centered 
at the following values of the gate voltage:E3ll3 
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where n is any integer. The peak splitting X+ — X^ 
grows linearly with Go . m 

In agreement with the experiments for a symmetric 
system, a = A = 0, the peaks at small conductance Gq 
are doubly degenerate (assuming Ec =^Ec, which is a 
good approximation for the experimentu). Even a small 
asymmetry, a <^ I, lifts this degeneracy. Indeed, the po- 
sitions of the peaks as Go ^ can be found from (H) 



with the electrostatic energy (||) as the full energy. As a 
result we find the two sequences of peaks 
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where again n is any integer. An asymmetry of the sys- 
tem caused by a non-zero A in Eq. (^) also leads to the 
lifting of the degeneracy. 

One can easily see that the peak positions given by 
Eq. show periodic beats: near certain values of the 
gate voltage X the neighboring peaks come very close 
together — they are separated by a distance of order a — 
while between those values of X the peaks are separated 
by 5X ^ 1. The period of these beats is a^^. In the re- 
gions where the distance between the neighboring peaks 
predicted by Eq. (Q) is small, an additional splitting due 
to the quantum charge fluctuations caused by finite inter- 
dot conductance Go should be taken into account. This 
additional peak splitting can be found in the same way 
as the splitting (^) in the symmetric case. For the case 
E(j = Ec, the result is 5X = (In 2/7r2)/iGo/e2. 

In the opposite case of strong coupling the properties of 
the system depend on the particular model of the junc- 
tion between the dots. For an electrostatically created 
constriction between the dots, a one-dimensioa.al (ID) 
model of the junction is the most appropriateEj In this 
case the conductance Go never exceeds 2e^//i, and the 
strong-tunneling case corresponds to a small reflection 
coefficient TZ = 1 — hGo/2e^ <C 1. We will concentrate 
on the asymmetric case, a > 0, assuming for simplicity 
A = 0, and derive the peak positions X* from Eq. (||). 
At fixed Af = Ni + N2 the electrostatic energy (^) can 
be rewritten as 



U^r{Nl) = Ec{Af-2Xf + 4Ec{Ni - 7)^ 



(8) 



where 7 = {Af + aX)/2. The second term on the right- 
hand side is expressed in terms of the number of parti- 
cles in the left dot. In the strong tunneling case Ni is no 
longer quantized, and at 7?. — > 1 its average assumes the 
value (A^i) — 7, thus minimizing the electrostatic energy. 
In this limit one easily finds Ej^{X) = Ec{Af ~ 2X)'^, 
and the peaks are equidistant, X* ~ {2n -f l)/4. 

At non-zero TZ the average number of particles in the 
left dot (iVi) is not precisely equal to 7, but oscillates 
near 7 with period A7 = 1. The corresponding small pe- 
riodic contribution to the ground state energy was found 
in Ref. |l^, where a single quantum dot connected to a 
large lead was considered. At temperatures exceeding the 
level spacings in both dots the tm) problems are equiv- 
alent, and one can use the resulted for the periodic cor- 
rection to the ground state energy. 
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where C = 0.5772 ... is Euler's constant, and 0(7) is 
defined as 0(7) = 00 + t^I- In contrast to Ref. |l^, 
we included here a shift (/iq in the phase 0(7) of the 
Coulomb blockade oscillations. Such a phase shift is al- 
ways present, for instance, due to an asymmetry of the 
scattering potential in the constriction connecting the 
dots—. In the case of a single dot connected to a large 
leadEl the presence of 4>q is irrelevant, as it can always 
be compensated by an appropriate shift of the gate volt- 
age. Similarly, in the case of a double-dot system the 
phase (pQ can be incorporated in the definition of X as 
some shift [see the definition of 7 in Eq. (||)], unless the 
system is completely symmetric, a — Q. 

One can now use the expression for the ground state 
energy (j^) to find the corrections to the equidistant peak 
positions caused by the weak scattering in the constric- 
tion. From Eq. (j5|) we find 



2n 
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+ (-1)" — --^7^1n — cos(20o + T^aX). 
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This result for the peak positions in the strong tunnel- 
ing regime, 7?. <C 1, is a generalization of the results in 
Refs. |l2|,^ to the case of an asymmetric system. The 
asymmetry gives rise to the cosine factor in Eq. (10). 



Similarly to the regime of weak tunneling, in the asym- 
metric case the distance between the peaks shows beats, 
with the period in X being a~^. 

As we mentioned, the presence of even a weak asym- 
metry of the system destroys the periodicity of the peak 
positions and thus complicates the comparison of the-epLi 
perimentally observed peak splitting with the theory£3'E£l 
One should note, however, that in both weak and strong 
tunneling cases, in the regions of the gate voltage X 
where the peak splitting assumes the smallest possible 
values, the distance between the neighboring peaks coin- 
cides with that predicted by the theoryll3'E3 for the sym- 
metric case. 

In the next section we calculate the heights and shapes 
of the conductance peaks, whose positions are given by 
Eqs. (^, (|^), and (p^), and compare the results with the 
available experiments. 



III. HEIGHTS AND SHAPES OF THE 
CONDUCTANCE PEAKS 

A. Weak tunneling between the dots 

We start our discussion of the heights and shapes of 
the conductance peaks with the case of weak tunneling 
between the dots, which means that the conductance of 
the constriction is small. Go <C /h. Nevertheless, we 
assume that the coupling to the leads is even weaker, 
Gi,Gr ^ Go. The results for the conductance depend 
on the symmetry of the system. In the symmetric case. 



one can find the conductance within the master-equatio, 
approach, identical to the one used in a single-dot case.l 
The resulting conductance has the form: 
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where (3 — AEq/T, and X'^ are the positions of the two 
adjacent peaks given by Eq. (^ with the same n. The 
two peaks are resolved only at sufficiently low tempera- 
tures, T < EciXl -XI) as shown explicitly in Fig. |. 
The derivation of Eq. (^l]) is outlined in Appendix 

In a symmetric device each state with an odd charge 
Af is doubly degenerate: the "odd" electron may be on 
either the left or right dot. In addition, at special val- 
ues of gate voltage X = X^. , states with charges A/" and 
JV+1 are degenerate. At these values of A, a charge can 
be transferred through the double-dot system via a se- 
quence of real states. As a result, the peak conductance 
is temperature-independent. 

A small asymmetry changes the situation qualitatively. 
As we saw in Sec. ||, the presence of a non-zero a or A 
in the electrostatic energy d) hfts the de generacy; i.e., 
the states with an "odd" electron on either the left or 
right dot now have different energies. If this difference 
of energies is larger than the temperature, one can no 
longer transfer charge through the double-dot system via 
real states only. For instance, at the value of the gate 
voltage A* = 1/(2 -I- a) the energies of the states with 
charge and with extra charge e on the left dot are 
equal, whereas the energy of the state with charge e on 
the right dot is A = 4:aEc/{2 + a). (Here we assume 
that all the asymmetry is due to a > 0, and A = 0, 
Ec — Ec-) At low temperatures the tunneling of an 
electron from the left dot to the right one involves an 
increase of the energy of the system by A > T, and so 
tunneling through the double-dot system is suppressed. 
Nevertheless, an electron can still escape from the left dot 
to the right lead via a virtual state in the right dot. Such 
a mechanism of tunneling in Coulomb blockade devices 
is known as cotunnelingH At temperatures exceeding the 
level spacing in.tlir.dot, the leading mechanism is inelas- 
tic cotunneling.EZllla In this case an electron tunnels from 
the left dot to the right one, and then another electron 
tunnels from the right dot to the right lead. After the co- 
tunneling process is completed, the right dot is returned 
to the state with no extra charge, but an electron-hole 
pair is created. As a result the phase space volume for 
such processes is proportional to T^, i.e., the conductance 
of the system is now suppressed at low temperatures. 

Let us find the height and shape of a conductance peak 
at low temperatures, T ^ A, for the example described. 
We first calculate the rate of cotunneling of an electron 
from the left dot to the right lead: 
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(12) 



where tpktsq/ ^ is the second-order matrix element for 
the transfer of an electron from state k in the left dot 
to state p in the right dot, and then the transfer of an- 
other electron from state q in the right dot to state s in 
the right lead; n^p q^g^ and ek(p.q,s) are the correspond- 
ing Fermi occupation numbers and energies, respectively. 
We also defined 

£ = C/(l, 0) - [/(0, 0) = 2(2 - a)Ec{Xl - X). (13) 

A straightforward calculation now yields 
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r(e) 3e' 
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e[l + {e/2TTTf] 
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t/T 



(14) 



Here we used the definition of the conductance Gq — 
/K)^ |tpfeP<5(efc)<5(ep), and a similar relation for Gr- 
One can now express the current through the system 

as 
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where wq and wi are the occupation probabilities of 
states with the charge of the left dot and e, respec- 
tively. Since the escape rate to the right electrode is 
strongly suppressed and is much smaller than the rate of 
tunneling to the left lead, the left dot is in equilibrium 
with the left lead and 
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Wo = 



l + g(e-ey)/T' 



Wl = 1 — Wo- 



(16) 



Here V is the bias applied to the leads. An expansion of 
the current (|l^) to linear order in V gives the conduc- 
tance 



{e/T)[l + {e/2nTr^ 
sinh(e/r) 



(17) 



The dependence e{X) is given by Eq. (|13|). 

As expected the height of the peak is suppressed at low 
temperatures as T^. The result ( |l7| ) can be applied to 
any of the peaks (|^) in the asymmetric system, provided 
that the appropriate values of A and e are found from 
the electrostatic energy (^. The cotunneling peaks cal- 
culated for realistic parametersQ are presented in Fig. |[ 



B. The limit of strong tunneling between the dots 



In section HI A we assumed that the coupling of the 



two dots is weak. This enabled us to apply the standard 
master equation technique in the case of a symmetric sys- 
tem, and to account for the lowest-order cotunneling pro- 
cess only in an asymmetric double-dot device. In this sec- 
tion we consider the limit of strong tunneling, where the 
techniques based on a simple perturbation theory are not 
applicable. We define the strong-tunneling limit as the 
case of perfect transmission through the channel between 
the dots, i.e., Gq — 2e'^/h. To treat this limit, we apply 
here a non-perturbative approach based on the bnsnpiaed 
picture of the ID transport through the channel.llSlij^E] 

We shall treat the double dot system as a single con- 
ductor of complicated shape. To find the conductance we 
will generalize the master-equation technique of Ref . ^ to 
account for the impedance of the charge propagation be- 
tween the dots due to the narrow constriction. We need 
to find the renormalization of the rates of tunneling from 
the leads into the double-dot conductor. The impedance 
of the charge redistribution within this conductor sup- 
presses the tunneling rates, not unlike the effectcj of the 
"electromagnetic environment" on transport through a 
single tunnel junction. 

To find the rate of tunneling through, e.g., the left tun- 
nel junction, we introduce a Hamiltonian that accounts 
for electron states in the left lead and left dot, as well as 
for the electron states participating in the re-distribution 
of the charge between the dots. The separation of the lat- 
ter group of states from the two others is possible at time 
scales shorter than the time of electron propagation from 
the tunnel junction to the other dot. In the case of a 
single mode constriction, this time is of the order of the 
inverse level spacing in the dot. Therefore our theory is 
limited to temperatures exceeding the level spacing. 

We assume that the constriction connecting the two 
dots is a single-mode channel with no reflection. In this 
case the set of electronic states responsible for the trans- 
port between the dots is one-|dipicnsional and can be pre- 
sented in a bosonized form.tSlla Thus the Hamiltonian 
can be written a,s H = Hq + He + Ht , 



He = EciriL - 2X)2 + AEc { m[u^{Qi) + ux(0)] - 
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Ht \ tkpa\apF + tlpalakF^j . 

kp 



dx 



-X 



(18) 

(19) 
(20) 
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First, in Hq, since we are considering transport of an elec- 
tron from the left lead into the left-hand dot, Ofc is the 
annihilation operator for electrons in the left lead and 
Up is the operator for electrons in the left dot; efc and 
Cp are the corresponding energies. The bosonized ID 
electron system is described by the displacements Ua-{x) 
and momentum densities Pcr{x) in two spin channels, 
which satisfy the commutation relation [ua-{x),Pa'{y)] — 
ihd{x — y)6a(T'', ^ and uq are the mass and density of 
ID electrons. Second, in Hq, the charging energy (^) 
is written in terms of the operator of the number of 
electrons tunneled through the left barrier and the charge 
eno[u|(0) -I- Mj.(0)] transferred from the left dot to the 
right one. Finally, in the tunnel Hamiltonian ( pO| ) the 
matrix elements tkp describe tunneling through the bar- 
rier. The transfer of each electron into the dot changes 
ul by one; to account for this, we use the operator F 



defined by the commutation relation 

[F,nL]=F. (21) 
The tunneling current through the junction is 

kp 

(22) 

where the average is performed with the density matrix 
of the system described by the Hamiltonian (|l8|)-(p0|). 
Assuming that the transmission coefficient of the tunnel 
barrier is small, we will calculate the tunneling current in 
lowest (second) order perturbation theory in tkp. Thus 
we can expand the density matrix up to the first order 
in tkp and find 



2e 



kp\ 



kp 



dt 



(at(0)a,(i))K(0)4(i))(i^t(o)F(t)) - (4 (i)a,(0))(a,(i)at(0))(i^(i)Ft(0))l • (23) 



In thermodynamic equilibrium the two contributions in 
II compensate each other. To find the effective conduc- 
tance Gl of the left junction, which is renormalized due 
to the slow charge redistribution between the dots, we 
now shift the chemical potential in the lead by eV and 
find Gl = diL/dV in the following formio 



Gl — ~r2^'^^ 



tK(t)dt 



-oc sinh [TTT{t — iS)/h] 



(24) 



Here G/ = (27re^/?i) J2 KfcpP'5(efc)'5(ep) is the unrenormal- 
ized conductance of the left barrier, and we have intro- 
duced the correlator 



K{t)^{F{t)F\0)). 



(25) 



In the derivation of Eq. (^J) we used the equality 
(F(t)Ft(O)) = (Ft(t)F(O)), which follows from the sym- 
metry of the Hamiltonian (^8|)~(|2^) with respect to the 
transformation — > —ul, F ^ F\ u„ —Ua, and 

In the absence of interaction, Ec — Ec — 0, the op- 
erators F and commute with the Hamiltonian, the 
correlator K{t) — 1, and the conductance is not renor- 
malized: Gl = Gl. We show below that the time depen- 
dence of the correlator K{t) is non-trivial if Ec, Ec > 0. 
Consequently the effective conductance Gl is renormal- 
ized and acquires a power-law temperature dependence at 
T<^Ec. 

To calculate K{t), we use a unitary transformation U 
which shifts the origin of the electron liquid displacement 
Ucr by a distance which depends on ul, 



U = exp 



1 + A 



-UL e 



(26) 



e 
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[p^{y) +Pi{y)]dy. 



(27) 



Upon the transformation ( |26| ) the Hamiltonian is simpli- 
fied, and operator F acquires a phase factor: 



U\Ho + Hc)U = Ho + EciuL - 2Xf 

+AEcnl[u^{Q) + uim\ 
.1 + A, 



tl^FU = F exp —I 



-e 



(28) 
(29) 



The correlation function ( p5[ ) now factorizes, K = 
KpKe. The factor Kp = (F(f)Ft(O)) is easily found: 



Kf 



-iEc{2nL + l-iX)t/n\ 



^-i4Ec(X'-X)t/h 
^4EciX'-X)/T i' 



(30) 



In the derivation of Eq. ( ^ ) we assumed that the gate 
voltage X is close to one of the peak positions X* = 
(2n -t- l)/4, see Eq. (p^), and that the temperature is 
much smaller than Ec so that only two states are in- 
volved. The calculation of Kq is also straightforward, 
as the Hamiltonian ( p^ ) is quadratic in the bosonic vari- 
ables, and the exponent in ( p9|) is linear in these variables: 



Ke{t) = { exp 



exp 



2 

(1 + A)2 



exp 



[e(o)-e(t)]e(o)> 



2ie'^Ec smh[nT{t - i5)/K\ 



One can now substitute K = KpKq into Eq. (^ 
find the renormalized conductance. 



(31) 



to 
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Gl = 



9i 
2 



4EciX - X*) 
T 



(32) 



where = (1 + A)^/4. The peak shape is given by the 
function Fr^{x) defined as 



F^{x) 



1 



r 1 + ? 



(33) 



cosh(a:/2) r(2 + ry) 

The tunneling into the double-dot system is suppressed 
at low temperatures, Gl oc T'^i^ . The origin of this sup- 
pression is Anderson's orthogonality catastrophe. The 
tunneling of an electron into the left dot results in a sig- 
nificant change of the ground state of the double-dot sys- 
tem, and the new ground state is orthogonal to the old 
one, leading to the suppression of tunneling. Indeed, af- 
ter the tunneling process has changed the charge of the 
left dot by e, charge q-\ = = e(l + A)/4 must be trans- 
ferred to the right dot in each spin channel to minimize 
the electrostatic energy. The orthogonality of the two 
ground states results in a power-law suppression of the 
tunneling density of states, G oc T''t_ where the exponent 
can be related to the charges a^ 77^ — '^Yl,ail<^ I ^)'^ ^ 
in agreement with Eq. (js^). 

The tunneling through the right barrier can be treated 
in the same manner, and the result for the renormalized 
conductance Gr can be found by replacement Gi Gr 
and ^ t;^ = (1 — A)^/4 in Eq. (^). After the renor- 
malized conductances Gl and Gn are found, we can use 
the master equation approach similar to the one outlined 
in Appendix M, and find the total conductance: 

G = (34) 

At T ^ the smallest of the two conductances, Gl and 
Gr, controls G, which means that the peak value of G 
is proportional to T"^ , with 77 = (1 -I- |A|)^/4. Depend- 
ing on the geometry of the system, the parameter A may 
vary from —1 to 1, and is in the symmetric case. There- 
fore the exponent of the temperature dependence rj varies 
from 1/4 in the symmetric case to 1 in the most asym- 
metric case. 



found that for a symmetric system in the weak tunnel- 
ing limit the conductance peak heights are temperature 
independent, whereas in the strong tunneling limit the 
peak heights are suppressed as T^/'*. For the asymmet- 
ric case we discovered the dependence of the peak 
conductance in the weak tunneling case, and with 
geometry-dependent exponent, \ < rj < 1, for strong 
tunneling. In this section we show that in the intermedi- 
ate regime the power-law temperature dependence of the 
peak conductance persists and consider the correspond- 
ing exponents. 

We start with the case of symmetric geometry. The 
weak tunneling result ( |lTl) was obtained in the frame- 
work of the master equation approach with the inter- 
dot tunneling rate calculated to first order in Gq. It 
is knownjlil however, that the higher-order terms of the 
perturbation theory give rise to a logarithmic renor- 
malization of the conductance Gq. This renormaliza- 
tion becomes important at temperatures T ^ Tk — 
Ec eyLvl-i-K^e^ / mGQf/\ Therefore at Gq < e^/h the 
result (O) is applicable only in the range of temperatures 
Tk-^T<^ Ec- 

A similar argument can be applied in the vicinity of 
the strong tunneling limit. Indeed, it was showrJij that at 
Go = 2 (e^ /ft.) (1 — 7?.) a weak refiection in the constriction 
is a relevant perturbation which grows at low tempera- 
tures and becomes strong at T < EcTZ- Thus the 



dependence of the peak conductance found in Sec. |III B 
holds only in the temperature range EcTi <C T <C Ec- 

To consider the effect of weak backscattering TZ I 
in the constriction on the low-temperature asymptotics, 
T <^ EcTi-, of the peak conductance, we complement 
the Hamiltonian (p^-(|20|) with a scattering term H' . In 
bosonic representation H' has the formlla 



H' 



D 



cos[27rrioUo-(0) — 0o] 



(35) 



C. Intermediate strength of tunneling between the 
dots 

In the previous two sections we considered the cases 
of weak and strong tunneling between the dots. We 



where the phase shift 0o is added to account for the pos- 
sibility of an asymmetric location of the scatterer with 
respect to the center of the constriction. One can then 
repeat most of the discussion of Sec. IIIB with the new 
Hamiltonian. Upon the unitary transformation (p6| ) the 
backscattering term takes the form 



27rnoUcr(0) - 0jc + -(1 + A)nL 



TT ^ 

a 

= - — \/^cos |7mo[uT(0) + Wi(0)] - (j^x + -^(1 + A)^^} cos{7mo[-ut(0) - mi (0)]}, 

TT y. 2 J 



(36) 
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where (px = (po + ^aX. 

Unhke other terms (^ ) of the Hamiltonian, the 
backscattering term (^) shows non-trivial dependence 
not only on the sum of the displacements u-f + ui, but 
also on their difference U| — mj^. At low temperatures 
T ^ Ec, one is only interested in the low-energy behav- 
ior of the system. In this regime the fluctuations of the 
charge of the dot no[u|(0) + mj (0)] are frozen due to the 
charging energy term in Eq. (Eq) and can be integrated 
out. The resulting backscattering term has the form 



U^H'U 



ie^EcD-R 



cos U)x - f (1-1- \)nL\ 



X cos{7mo[MT(*-') ~ "i(0)]}; 



(37) 



cf. Ref. ^ Since the operators F and F'^ do not com- 
mute with the backscattering term (|37|), the latter can 
affect the Kp component of the correlator K{t) and the 
conductance of the left barrier (|4|). To find the ef- 
fect of the backscattering on Kpit), we first discuss the 
influence of the operator (|37| ) on the dynamics of the 
spin field U| — mj^. One can easily shew that the op- 
erator ( ^ ) is a relevant perturbation,c2l i.e., the ampli- 
tude of the cosine term grows at low energies. Thus at 
T — > the spin field fluctuations are frozen at the value 
no[u|(0) — uxC^)] = or 1 for the positive and negative 
values of cos ^(px — f (1 + A)ni] , respectively. 

When an electron tunnels into the double-dot system 
through the left barrier, the value of ul changes from 
to 1. Thus the prefactor in Eq. ( ^7| ) is proportional 
to either coscpx or cos [0x — f (1 + A)ni,]. If the two 
cosines have the same sign, the increase of n^, described 
by the operator does not affect the long-time dy- 
namics of the spin mode, which remains pinned at the 
origin with the same value of M|(0) — Ui(0). In this 
case the time dependence of Kpit) is not affected by 
the backscattering and the conductance Gl is still given 
by Eq. (^2|), with a different prefactor which we do not 
calculate here. On the other hand, if the signs of cos (px 
and cos [(px — f(l + A)nL] are different, the change of 
TT-L shifts the boundary condition for the spin mode from 
no[w|(0) — it|(0)] = to 1. An abrupt change of the 
boundary condition creates a disturbance in a ID bosonic 
field that decays slowly, giving rise to power-law time 
dependences of electronic Green's functionsBj Thus the 
correlator if J? (t) acquires an additional time-dependent 
factoiia^y KJt) cx nT / i smh[T:T {t ~ iS)/h]. Accord- 
ing to Eq. (|24|), such a modification of K{t) not only 
changes the prefactor in Eq. (p2|), but also replaces the 
exponent 77^ by r;^ -|- 1. Thus depending on the values 
of (px and A, the backscattering either does not affect 
the temperature dependence of the renormalized conduc- 
tance Gl oc t(^+'^) or replaces it with a stronger one, 
Gl oc ri+(i+A)V4. 

The latter result can be easily interpreted in terms of 
the orthogonality catastrophe. Indeed, as we saw, the 
tunneling of an electron into the left dot leads to the 



transfer of charge + qi = e{l + X) /2 through the con- 
striction. On the other hand, if no[u|(0) — M|(0)] changes 
from to 1, the transferred spin is {q^ — qi)/2e — i. Thus 
we can easily evaluate the charge transferred in each of 
the spin channels. 



= e ± 



1 l + \ 



(38) 



The suppression of the tunneling density of states is de- 
scribed by|-tlie power law v cx e''^-, where the exponent 
is given byn3 77^ = 2J2ail<^/^)'^- From Eq. ( |3^ ) we now 
find 



'7l = 1 



(39) 



which results in the power-law suppression of the con- 
ductance, Gl oc T^+ii+^r/\ 

To find the temperature dependence of the peaks in 
the conductance through the double-dot system, one has 
to find not only Gl, but also G^. It is clear that when a 
tunneling process through the whole double-dot system 
is completed, exactly one electron is transferred through 
the constriction. Thus we conclude that the total transfer 
of charge after one electron has tunneled into the left dot 
and another one escaped from the right dot is Aq — e. 
Since the charge transferred through the constriction at 
the first step was e(l + X)/2, the tunneling of an electron 
from the right dot must be accompanied by the transfer of 
charge e(l — A)/2. We saw above that unlike charge, the 
spin is transferred in quantized portions As = ^. Since 
the total transferred spin is ^, we conclude that exactly 
one of the two tunneling events involved the transfer of 
spin. Therefore in the cases when the temperature depen- 
dence of Gl is given by T^-^+^l^/^ and ri+(i+^)V4^ the 
conductance Gr behaves as T^+i^-^f/'i ^nd T^^->^f/i^ 
respectively. Finally, since the total conductance is given 
by the smaller of and Gr, we find 



G< 



^^^CL_AL^ if cos(/,xcos[(/)x - f(l + A)nt] > 0, 
r^, if cos(/)jf cosUx — f (1 + A)n/^1 < 0. 



(40) 



To determine which option applies to a particular peak, 
one needs a detailed knowledge of the microscopic struc- 
ture of the double-dot system. It is clear, however, that 
for nearly symmetric geometries the parameter A, which 
is determined only by the electrostatics of the system, 
should be small: A ^ 1. In this case we predict the tem- 
perature dependence G oc T^/'* for all peaks, independent 
of the microscopic structure of the double-dot. 

The temperature dependence of peak heights has been, 
investigated experimentally by N. van der Vaart et ai.E£l 
In the regime of weak reflection, the power-law fit to the 
data in the temperature interval lOOmK < T < IK gave 
77 = 0.8-1.2, which is slightly less than our result r] = 1.25 
for the symmetric geometry. 
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The result (gO|) shows that the presence of even weak 
backscattering in the constriction gives rise to a large cor- 
rection Ary ~ 1 to the exponent in the power-law temper- 
ature dependence G cx T"^ of the peak conductance. In 
the derivation of Eq. ( ^0| ) we assumed that the backscat- 
tering is weak, 7?. <C 1. As the backscattering grows, 
it further affects the temperature dependence. Indeed, 
so far we assumed that the presence of the backscatter- 
ing only creates a boundary condition for the spin mode, 
M| — Uj, and does not affecjt-^he charge mode, M| + u^. 
However, from the studieai3li3 of a single dot connected 
to a large lead it is known that the backscattering does af- 
fect the charge transferred through the constriction. One 
can attempt to generalize the result ( |40| ) to the case of ar- 
bitrary TZ by introducing the value Qt of the charge trans- 
ferred through the constriction after an electron tunnels 
into the left dot. It is clear from the derivation of Eq. (|39| ) 
that the second term there is actually Qt/e, i.e.. 



(41) 



At 7?. ^ the correction to the electrostatic value 
Qt = e(l 4- A)/2 of the transferred charge is small,t3 
AQt ~ TZln^^ which justifies the approximation (^9|). 

It is interesting to apply Eq. ( ^ll) to the weak-tunneling 
limit 7?, ^ 1, considered in Sec. Ill A. In the asymmetric 



case we demonstrated that the temperature dependence 
of the conductance is given by Eq. ([l7|). In the derivation 
we assumed that because of the high barrier separating 
the dots, there is no transfer of electrons between the 
dots after an electron tunnels into the left dot. This 
means that both Qt and As vanish, r/^ = 0, and the 
conductance Gl is not suppressed at T ^ 0. On the 
other hand, when an electron escapes to the right lead, 
it must go through the constriction, leading to Qt = e 
and As = \- As a result, the relation for 77^ similar to 
Eq. (^l|) will give ijj^ = 2, which leads to the quadratic 
temperature dependence ( p^ ) of the linear conductance. 

In our approach the shape of the peak is obtained 
from Eq. ^ by substitution K(t) cx {irT/i sinh[7rT(i - 
i5)/h]}^. As a result the shape of the peak is always 
uniquely related to its temperature dependence. 



G cx T^F, 



(42) 



where F,^ is defined by Eq. (33), and e is proportional 



to the deviation of the gate voltage from the peak cen- 
ter. One can easily check that the peak shape (|l^ ) in the 
weak-tunneling limit does coincide with F2{s/T). 



and thus affect the Coulomb blockade phenomenon that 
develops at a sufficiently low temperature, T < Ec- The 
positions of the peaks in linear conductance, G{Vg), de- 
pend on the value of Go, and in the limit Go — > 2e'^/h the 
peaks become equidistant (Section ||). A striking result, 
however, is that the height and shape of the peaks also 
evolve significantly with Go, and remain non-trivial even 
in the limit of a refiectionless constriction. Go 2e^/h. 
We have demonstrated that at any Go, except a special 
case of small Go in a symmetric two-dot system (Sec- 
tion [II A), the peak conductance is a power-law func- 



tion of temperature T. The exponent of the power law 
depends on the charge re-distribution between the dots 
that accompanies the electron transport through the two 
dots. Both this exponent and the explicit peak shapes 
depe nd o n the dots' geometry, as well as Go (Sections 
IIIB and III C ). The suppression of conductance at low 
temperature and bias can be understood in quite general 
terms as an Anderson orthogonality catastrophe caused 
by the re-distribution of charge (see Sections IIIB, 111 C ) , 
and the same exponents should describe the bias depen- 
dence of the differential conductance at T = 0. 

In deriving our results, we assumed that the incoming 
electron dwells in a dot for a long time td ^ fi/T be- 
fore reaching the constriction that connects the dots. In 
a generic situation of a dot lacking a special symmetry, 
the electron bounces off the walls many times, before it 
gets to the constriction. The dwelling time is determined 
by the level spacing, td ^ h/6E. Therefore, the results 
we presented in Section are valid in the temperature 
interval SE < T < Ec- For typical parameters, this in- 
terval allows one to vary the temperature by at least one 
decade. 

The effect of quantum charge fiuctuations on the 
ground state energy has been recently demonstrated 
experimentally.QQ The data of Waugh et alp is in a quan- 
titative agreement with the present theoryJlj The tem- 
perature dependence of the peak conductance, which is 
related to the dynamics of the charge re-distribution, was 
studied in a very recent experiment by N. van der Vaart 
et aZHa The temperature dependence exponent found ex- 
perimentally in the regime of weak reflection, rj = 0.8-1.2, 
is somewhat smaller than the theoretical value 77 — 1.25 
we find for this case. 
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APPENDIX A: DERIVATION OF THE 
ELECTROSTATIC ENERGY EQ. (§) 

In this section we find the electrostatic energy of the 
double dot structure in terms of the capacitances of the 
individual dots and the gate voltage. To describe the 
electrostatics of the physical structure shown schemati- 
cally in Fig. ^ we introduce the circuit diagram in Fig. ^. 
The electrostatics is determined by the gate voltage Vg 
and 5 capacitances — C'l and C2 are the capacitances of 
the dots to the gate, C3 is the capacitance between the 
dots, and C4 and C5 are the capacitances of the dots to 
everything else. 

In terms of the charge on each capacitor, the electro- 
static energy is 

(Al) 

The number of electrons on each dot is given by the sum 
of the charges on the appropriate three capacitors: 



Ec = 



eNi = qi+ q3 + (74, 
-eN2 = (?2 - 93 + 95- 



(A2) 
(A3) 



We now must minimize the energy Eq. (Al) at fixed val- 
ues of Vg, Ni, and N2 and evaluate the energy at this 
minimum. The result has the form (4|) [up to an irrele- 
vant constant], and we now give explicit expressions for 
the parameters in this equation. 

First, the energy involved in changing the total charge 
on the double-dot system is given simply by the total 
capacitance of the double-dot to the external world. In- 
troducing the external capacitance 



we find 



Cext = Ci + C2 + C4 + C5, 



Ec 



2Ce 



(A4) 



(A5) 



The coupling of the total charge to the gate is given by 
the capacitance to the gate. 



X = 



-{Cl+C2)Vg 

2e 



(A6) 



Turning to asymmetric structures, we find that the frac- 
tional asymmetry of the capacitances determines the pa- 
rameter A, 



A — (C2 + C5 — Ci — C4)/ Coxt- 



(A7) 



In terms of this asymmetry parameter, we find that the 
charging energy for transfer from one dot to the other is 



2[Cext(l-A2)+4C3] 



(A8) 



and that the coupling of this excitation to the gate is 
given by 



a = 2 X 



Ci — C2 
C1+C2 

This completely specifies the electrostatic problem. 



(A9) 



APPENDIX B: MASTER EQUATION 
TECHNIQUE FOR THE CONDUCTANCE PEAKS 



In this Appendix we derive the expression ( |11| ) for 
the conductance peaks in the symmetric case. We re- 
strict ourselves to the case of temperatures which are 
much smaller than Ec, but can be of the order of 
Ec{X*^ — Xl_). For simplicity we assume that the gate 
voltage X is close to | , so that only the pair of peaks cen- 
tered at and XI given by Eq. (||) with n = should 
be considered. In this regime only the states with charges 
0, e, and 2e should be taken into account. Clearly, one 
has four states, which can be denoted as 0, I, r, and 2, 
where / and r describe the two states with charge e on 
the left and right dot respectively. 

We start with introducing the probabilities of occupa- 
tion of the four states, which satisfy the obvious condi- 
tion: wo + wi+'Wr + W2 — 1. The rate of transitions from 
state to state I, which are caused by tunneling of an 
electron through the left barrier, is given by 

1 27r ^ — A . .0 , >^ 
= Wo— } \tpk\ nk(l - Up) 

kp 

xS{ek + eV + Eo~ep~Ei) 
= ^wofiEi-Eo-eV). (Bl) 

Here tpk is the matrix element of tunneling from the state 
k in the left lead to the state p in the left dot, 
and efc(p) are the corresponding Fermi occupation num- 
bers and energies, Eq{X) and Ei{X) are the values of 
the ground state energy of the double-dot system with 
charge and 1, the function f{x) is defined as 



(B2) 



In deriving Eq. (Bl) we assumed the bias V applied to 
the left lead. In a similar way one can find the rates of 
all other transitions. 

In a stationary state the time derivatives of the oc- 
cupation probabilities of all the four charge states of the 
system vanish. This gives us three independent equations 
wi = Wr — W2 = 0, which can be written in the following 
form: 
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Gi[wof{ei - eV) - wif{~ei + eV)] + Go/(0)K " wi] + G,K/(-e2) - wif{e2)] = 0, (3) 
G.K/(ei) - Wrfi-si)] + Go/(0)H - Wr] + Gi[w2f{-e2 + eV) ~ Wrf{s2 - eV)] = 0, (4) 
G,K/(£2 - eV) - W2f{-e2 + eV)] + Gr[wif{e2) - W2f{-e2)] = 0. (5) 

Here we have introduced 

ei = Ei{X) - Eo{X) = 4Ec{X*_ - X), 
62 = E2iX) - EiiX) = AEciXl - X). 

The current / can be also expressed in terms of occupation probabilities Wi. In a stationary state, currents through all 
the junctions are equal. Considering the current through the link between the dots, we can express / in the following 
form: 

I=^f{0){wi-Wr). (6) 

e 

The equations (|3|)-(|^) must hold at any bias. Since we are interested in linear regime eV ^ T, we can differentiate 
equations @~(^ and replace them by the system of equations for the derivatives w'q, w'l, w^, and w'2 of the occupation 
probabilities over bias: 

^)+Ro(^-^)+R.r(^~^) =0, (7) 

TJ \Wr WiJ \W2 Wl J 

f w'l w'\ f Wn w' e\ 

f i?o M - — + i?2/ M - ^ - 7^ = 0, 8 

\Wl WrJ \W2 Wr 1 J 

i?.f^-^ + |;)+i?2.f^-^) ^0, (9) 

\Wr W2 1 I \ Wl W2 I 



Ri 








Wl 


Rr 


\Wo 


w'r 

Wr 



cf. Ref. 27 Here Wt are the equilibrium occupation probabilities, and we introduced the equilibrium rates: 



Ri 


= woGif{ei) 


- wiGif{- 


ei)^Gif{ei)Z-\ 


Rr 


= woGrf{ei) 


= WrGrf{- 


-ei) = Grf{ei)Z-\ 


R{) 


- uvGo/(0) = 


= i«;Go/(0) 


^TGoe-''/^Z-\ 


R2I 


= WrGlf{£2) 


^W2Glf{- 


e2)^Gif{e2)e~'''^Z-\ 


R2r 


= WiGrf{£2) 


= W2Grf{- 


'e2) = Grf{e2)e-'''^Z-^ 



(10) 

(11) 

(12) 
(13) 
(14) 



with the equilibrium partition function Z = \ + 
2exp(— £i/r) + exp[— (ei +e2)/T]. Since the sum of oc- 
cupation probabilities always equals one, there arc only 
three independent variables in the system 

Using Eq. (0), we can also express the conductance 



G = i?o ^ - ^ (15) 

\Wl WrJ 

in terms of the solution of the system (0)-(^). In the 
limit Go ^ Gi,Gr, we find: 

Q _ }_ f RlRr _|_ R2lR2r \ ^-^^^ 
T \Ri + Rr R2I + R2r J 

Substitution of the rates (O)-(ra into (111) yields the 



formula (11) 
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FIG. 1. Schematic view of the double quantum dot system. The dots are formed by applying negative voltage to the gates 
(shaded); the solid line shows the boundary of the 2D electron gas (2DEG). Vi and Vr create tunnel barriers between the dots 
and the leads while Vb controls the transmission coefficient through the constriction connecting the dots. 



FIG. 2. The evolution of the split peaks with temperature described by Eq. (^l|). The reference point for the gate voltage 
is chosen to be (X* -|- X^)/2, and the gate voltage is plotted in units of (XjJL — XI) /2. The peak splitting is observable at a 
sufficiently low temperature, (3 = AEc/T > 2. 



FIG. 3. The conductance as a function of dimensionless gate voltage X in the asymmetric weak-coupling case. Note the 
correlation between the modulation of the peak height and the separation of adjacent peaks: when the peaks are high the 
splitting is small, while when the peaks are small they are well separated. The parameters used in Eq. (nTl) to produce this 
plot are a = 0.155 and T/Ec = 0.07. 



FIG. 4. Equivalent electrostatic circuit for the double-dot device of Fig. hi in equilibrium. 
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Fig. 3 [Matveev, Glazman, and Baranger] 




Fig. 4 [Matveev, Glazman, and Baranger] 



